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Abstract 

We compute bounce solutions describing false vacuum decay in a model 
in four dimensions with quantum back-reaction. The back-reaction of the 
quantum fluctuations on the bounce profiles is computed in the one-loop and 
Hartree approximations. This is to be compared with the usual semiclassical 
approach where one computes the profile from the classical action and de- 
termines the one-loop correction from this profile. The computation of the 
fluctuation determinant is performed using a theorem on functional deter- 
minants, in addition we here need the Green' s function of the fluctuation 
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operator in oder to compute the quantum back-reaction. As we are able to 
separate from the determinant and from the Green' s function the leading 
perturbative orders, we can regularize and renormalize analytically, in anal- 
ogy of standard perturbation theory. The iteration towards self-consistent 
solutions is found to converge for some range of the parameters. Within this 
range the corrections to the semiclassical action are at most a few percent, 
the corrections to the transition rate can amount to several orders of magni- 
tude. The strongest deviations happen for large couplings, as to be expected. 
Beyond some limit, there are no self-consistent bounce solutions. 
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1 Introduction 



False vacuum decay |2l El is one of the basic mechanisms which are of- 
ten invoked in the construction of cosmological models or scenarios. False 
vacuum decay may initiate inflation, or it may happen after inflation, if the 
universe gets trapped in one of the local minima of the Higgs potential in 
Grand Unified Theories. The false vacuum may decay by spinodal decom- 
position j3J in] if the minimum becomes a maximum by a change of the 
effective potential, e.g., with decreasing temperature or cosmological expan- 
sion, or triggered by another field (inflaton) [7j . If this does not happen it 
may decay by an over-the-barrier transition at finite temperature, the process 
of bubble nucleation, or it may tunnel to the true vacuum, a process that is 
described semiclassically by bounce solutions of the associated Euclidean the- 
ory. If gravity is included, such a solution is the Coleman-de Luccia bounce 
[HI E] • We will concentrate here on bounce solutions in a theory with one 
single scalar field and an asymmetric double-well potential. The methods 
described here may be transferred to more realistic models. In particular, 
they may also be used for the multifield case with coupled equations for the 
classical and fluctuations fields. 

The most common approach to computing the tunneling rate is the semi- 
classical one. The bounce profile is determined by minimizing the classical 
Euclidean action while in the calculation of the transition rate one includes 
the one-loop prefactor, i.e., the fluctuation determinant computed around the 
classical profile. For false vacuum decay in 3 + 1 dimensions this computa- 
tion to one-loop accuracy has been presented in Ref. [10], based on methods 
developed in ^T]. Recently similar computations have been presented in Ref. 
fr2\ using ( function regularization. A different technique for computing the 
fluctuation determinant has been used in [T31 CH US] . 

Once the quantum corrections become important one may ask whether 
these quantum fluctuations react back upon the bounce profile. Two ways 
of including the quantum backreaction for such semiclassical systems have 
been discussed in the literature: the one-loop backreaction and the Hartree 
backreaction. The backreaction to one-loop order consists in determining 
the classical field in such a way that it is an extremum of the one-loop 
effective action, the sum of classical and one-loop action. This has been 
discussed by Surig [12] in the context of bubble nucleation in a SU (2) gauge 
theory (a multifield case). We will call this the one-loop backreaction. In 
the Hartree approximation one takes into account the backreaction not only 
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upon the classical field, but also upon the quantum fluctuations. We will 
call this the Hartree backreaction. This approach has been applied recently, 
using quite different mathematical and numerical approaches, by Bergner and 
Bettencourt |12j and by the present authors [18j to the case of false vacuum 
decay in two space-time dimensions. Here we extend the latter scheme to 
3 + 1 space-time dimensions, i.e. four Euclidean dimensions for the bounce. 
We also present results for the one-loop backreaction. 

The methods used for the one-loop and the Hartree approaches are closely 
related. In both cases one has to compute not only the fluctuation determi- 
nant but also the Green function Q{x,x') of the fluctuation operator. This 
is done in analogy to the two-dimensional case ^H] • For the case of a bounce 
in four Euclidean dimensions renormalization becomes more involved than 
the one in two dimensions, by the occurrence of logarithmic, quadratic and 
quartic divergences, both in the Green's function ai x = x' and in the fluctua- 
tion determinant. So regularization and renormalization have to be discussed 
more extensively. Besides the regularization technique we will have to dis- 
cuss the renormalization conditions. One way of renormalizing is to use the 
MS prescription on the basis of dimensional regularization. This implies a 
rather transparent choice of counterterms and it would be appropriate if the 
bounce were computed for a theory, in which the parameters are determined 
from experiment, using the same prescription. This prescription has been 
used in Refs. [TUl CHI On the other hand, in our toy world, one may 
wish to compare the results of the various approaches for situations were the 
effective potentials have similar "physical" properties, e.g. the same position 
of the local minima and the same energy differences between the two vacua. 
This approach is intuitively more appealing, while the determination of the 
counterterms gets somewhat less transparent. These renormalization condi- 
tions were used in Ref . jTT] . In the latter publication the authors have used 
Pauli-Villars regularization. This could have been done here as well; as the 
divergent parts are separated in a Lorentz-invariant way from the flnite parts 
which are computed numerically, the method of regularization can be chosen 
freely, as in usual perturbative calculations. 

The plan of the paper is as follows: In section |21 we present the basic 
equations of the theoretical model and for the tunneling rate. In section 
El we introduce the formalism used to compute the quantum backreaction 
in the Hartree approximation, using the 2PPI expansion for the effective 
action. In section |3] we present the formalism for computing the quantum 
backreaction to one-loop accuracy, as the lowest order of the IP/ action. We 
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then turn to the analytical and numerical methods. In sections and El we 
introduce the mathematical formalism used for the numerical computation 
of the Green's function and of the fluctuation determinant, respectively. In 
section [7| we discuss the particularities associated with the unstable mode in 
the n = and with the translation mode in the n = 1 partial waves, and the 
numerical techniques used to overcome these problems. The renormalization 
of the dynamics is discussed in section |H1 the technical details are deferred 
to Appendices El and El The numerical results are presented and discussed 
in section We summarize our approach and the results in section ITUl 

2 Basic relations 

We consider a scalar field theory in 3 + 1 dimensions, with the Lagrange 
density 



and is represented in Fig. ^ (solid line). It displays two minima, one at 
$ = 0, representing the false vacuum, and the other one at $ = (ptv > 
corresponding to the true vacuum. 

The classical Euclidean action of this theory is given by 



(2.1) 



The potential U ($) is given by 
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(2.2) 




(2.3) 



and the bounce which minimizes this action satisfies 




(2.4) 



or, using its spherical symmetry. 



(f) 3d(j) 



+ U\<P) = . 



(2.5) 



r dr 
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The boundary conditions are 

d(j) 

■^|r=0 = 0, 0r-^oo = (Ptv ■ (2.6) 

The one-loop correction to the classical action is given by 

1 -d^ + U"(^(x)) 1 

Si-i = 2 l^det ^ = 2 l^^^M ' (2-7) 

where is the Laplace operator in four dimensions, m is the mass in the 
false vacuum, 

= f/"(0) . (2.8) 

and where the prime indicates that the translation zero mode is removed 
and that one replaces the imaginary frequency of the unstable mode by its 
absolute value. 

The transition rate from the false to the true vacuum is given, in the 
semiclassical approximation without backreaction, by 

r^-^°°'^ = (1^) ' V~'/' exp(-5e,) = ' exp(-5,i - S,.,) , (2.9) 

where 

m = det' Y ^ . (2.10) 

We recall (see e.g. 19 ) that the prefactor arises from the quantization 
of the collective coordinate associated with the translation zero mode. The 
presence of a zero mode rjo in this approximation is demonstrated by taking 
the gradient of the classical equation of motion: 

V, [-9^0 + [/'(0)] = [-d^ + f/"(0)] V.0 = . (2.11) 

Its normalization, defined by rjoi = MoViCp, and the condition that ?7oi is 
normalized to unity, is given by 

J^~^ = J d^x{y4f = T/2, (2.12) 

where there is no summation over i and where T is the kinetic part of 
the action. Quantization of the collective coordinates yields the prefactor 
(A/'o"V2vr)2 = (T/47r)2. Furthermore, one finds in four dimensions, using a 
scaling argument, that V = J d'^xU{(f)) = —T/2. We then obtain T = 2S'ci 
and the prefactor becomes (T/47r)^ = (5'ci/27r)^, as written down in Eq. 

(ESD- 
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3 The bounce in the Hartree approximation 



The Hartree approximation can be derived in various ways, in the most 
intuitive approach from a variational principle using a wave function which 
is a direct product. Within the so-called 2PPI formalism [201 El 1221 
lowest quantum approximation. In terms of Feynman graphs it represents a 
resummation of daisy and super-daisy diagrams, like in the large- iV effective 
action. The effective action in this formalism is given by 

up to renormalization counterterms discussed in section |H| Here A is a local 
insertion into the propagator which has the form 

g-\x) = -d^ + M\x) , (3.2) 

with the definition 

M^ = m^- 6r/0 + + -AA = t/"(0) + -AA . (3.3) 
A itself is defined by the equation 

With this definition Eq. ()3.3p becomes a self-consistent equation, the gap 
equation. Finally r^^^^ is the sum of all two-particle point-irreducible graphs, 
in which all internal propagators have the effective masses Ai^. A graph is 
two particle point reducible (2PPR) if it falls apart when two lines meeting 
at a point are cut. To lowest order in a loop expansion r^^^^ is given by a 
simple loop, i.e. 

r^^^^ = ^lndet- ; , (3.5) 

and this is equivalwent to the Hartree approximation. We will discuss the 
modifications arising from the zero and unstable modes later. In this lowest 
approximation, A is given by 

^^^^ " '^ SM^ix) ^ Lg2^^2 l^ >= ^(^'^) = •^(^) • (3-6) 
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Here the Green's function Q is defined by 



(3.7) 



and where we have introduced the fluctuation integral 



^{x) = Q{x, x) . 



(3.8) 



In taking variational derivatives of the effective action we have to consider 
A as a functional of Ai^ and (j), i.e., in the last term of Eq. ()3.1|1 we have to 
replace 



see Eq. (Q- 

Taking the variational derivative of the effective action with respect to 
A4^ leads back to Eq. ()3.3p . The partial derivative with respect to the field 
(f) leads to the equation for the bounce profile 



It is not a simple differential or integro-differential equation as jF(x) is a 
nonlinear functional of 0. 

Using rotational symmetry we obtain explicitly 



The backreaction of the quantum modes upon themselves is contained in 
J-'{x), or, equivalently, in J\4^{x). For a spherical background field, as we 
have it here, these functions are themselves spherically symmetric. 

Of course J-'{x) as the limit x' x of the Green's function Q{x,x') is 
ill-defined, as is the fluctuation determinant in Eq. ()3.5j) . This problem will 
be dicussed in section |H] 



m' + 6r]<p) = ^{M'-U"{<P)) , (3.9) 




(3.10) 
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4 The bounce with one-loop backreaction 

The one-loop backreaction can be derived using the IPI formahsm and rep- 
resents in fact a simphfication in comparison with the Hartree backreaction. 
The action is given by 

5'eff,l-l[>^^ 0] = SM + r^^^[0] , (4.i) 

up to renormahzation countcrtcrms. Here r^^-'^[0] is the sum of all one- 
particle irreducible graphs with the propagators 

g'^^x) = -d^ + M\x) , (4.2) 

where here Ai^ is given by 

M^ = m^ - 6r](p + ^A0^ = U"{(j)) . (4.3) 

To the lowest order is given by 

r- = i,„det'4±4. (4.4) 

Taking the partial derivative with respect to the field leads to the equation 
for the bounce profile 

+ V\<^{x)) + I \\<\){x) - 2r]] T{x) = . (4.5) 

Zi 

The last term in this equation arises by 

sm - TOW ^ 2^^^)^^^^^^) - ' ^'-'^ 

where again J-'{x) = Q{x.x). The similarities with the Hartree backreaction 
are obvious, we will need analogous numerical techniques, the main differ- 
ence are in the effective mass, which in the one-loop formalism is simply a 
function of 0, whereas in the Hartree approximation it is determined in a 
self-consistent way. 
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5 Computation of the Green's Function 



In order to include the backreaction of the quantum fluctuations on the 
bounce in the Hartree approximation we need the Green's function Q{x,x') 
of the (new) fluctuation operator. In fact the Green's function is usually 
discussed in a more general form, as a function of energy. We will need to 
take this more general approach in order to be able to discuss the translation 
mode and in order to discuss the determinant theorem in the next section. 
Such a concept corresponds to introducing an additional fifth dimension. We 
will choose it spacelike, thus introducing an additional Euclidean time. As </> 
still lives in four Euclidean dimensions, we have translation invariance in the 
new time direction. So we can introduce the Fourier transform ^(x,x', z/^), 
where u is the Euclidean frequency. The Green's function now satisfies 

[-9^ + m^ + V{r) + iy^]g{x, x' , u^) = 6\x - x') , (5.1) 

with 

V{r) = M\(l)) -m^ . (5.2) 

Here A1^(0) is given either by Eq. ()3.3|) or by Eq. ()4.3|) in the Hartree and in 
the one-loop formalism, respectively. In the one-loop formalism A^^(O) = m^, 
so with the boundary conditions (j2.6j) y(r) ^ as r ^ oo. In the Hartree 
formalism is, via the Green function, a functional of 0. Here should 
be equal to for a constant field 0(x) = 0. This has to be imposed as 
a renormalization condition. Alternatively, in the definition of V{r) one 
replaces the bare mass with the renormalized mass in the false vacuum. 

The Green's function can be expressed by the eigenfunctions of the fluc- 
tuation operator. We denote them by ipa{x), they satisfy 

[-d^ + m^ + V{r)]Mx) = ujIM^) ■ (5.3) 

In terms of these eigenfunctions the Green's function can be written as 

I 2\ V^a(2;)'?/'o(x') , , 

^x,x,z/ => . 5.4 

a " 

We may, furthermore, decompose the Hilbert space into angular momentum 
subspaces, introducing eigenfunctions Ynim{S^^)Rkij')^ where YnimiS^z) are the 
spherical functions on the 3-sphere (see e.g. the Appendix of Ref. ISHj) and 
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where the radial wave functions -R„Q,(r) are eigenfunctions of the partial wave 
fluctuation operator: 



3 d , nin + 2) , ^ 



-— + -^—y — - + + V{r) 



dr^ r dr r 



Rna{r) = ul^Rnair) . (5.5) 



The index a labels the radial excitations, the spectrum is continuous, but 
may include some discrete states, e.g., the unstable and translation modes. 
In terms of these eigenfunctions the Green's function takes the form 

nlm a "'"^ 

While these expressions are very suitable for discussions on the formal level, 
they are not very suitable for numerical computation. In particular, if one 
uses these expressions for the numerical computation, it becomes necessary 
to discretize the continuum states by introducing a finite spatial boundary. 

There is a well known alternative way of expressing Green's functions. 
Consider first the free Green's function obtained for V{r) = 0. It can be 
written as 

/d'^k ^ik-ix-x') 
/o / 2 I ^ — 1 ' (^•'') 

and this may be expanded as 

^4^ ^k-(x-x') ^ ^2 

(27r)4 p + m2 + z/2 ~ 47r2 kR 

^ ^( , ^\n\( ^/n+l(/^r<) ir„+i(/«r>) 
= ir^lX^^Wrk^^^X) 

n=0 ^ ^ 

^ oo 

= TT^ y^l'^ + l)Cn(cosx)xi°i(r<, fi;)/ci°|(r>, k) . 

n=0 

With r = |x| and r' = \x'\ we have r<; = min{r, r'}, r> = max{r, r'}, k is 
defined as k = ^rrP- + z/^. R and x ^i-re defined by 

Q I / 1 2 2 /2 / / \ 

i? =|x — x| =r +r — 2rrcosx, (5-8) 

i.e., X is the angle between the directions fis and fig of x and a;'. The 
functions C\ are Gegenbauer polynomials, see section of Ref. jSlj. The 
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expansion of Ki{k,R) in terms of products of Ii{nr^) and Ki{nr^) is the 
Gegenbauer expansion, given in section 7.61, Eq. (3) of Ref . |24j . For the 
case X = x' one has C^(cosx) = = n + 1. For convenience we have 

introduced the functions X„ and Kn as 



J„(r, k) 
/C„(r, k) 



'n+l 



n+l 



'Kr)/r , 



They satisfy 



3 (i nin + 2) , 



i3„(r, k) = , 



(5.9) 



(5.10) 



(ir^ r dr 

where Bn stands for X„ or /C„. Their Wronskian is given by 

/C„(r, K)dln{r, k) / dr - In{K,r)dlCri{r, k) / dr = 1/r^ . 
We expand the exact Green's function in an analogous way with the ansatz 

^ oo 

G{x, x\ u') = ^,Y.^n + l)C^(cosx)/-(r<, «:)/+(r>, k) . (5.11) 



n=0 



The functions /^(r, n) satisfy the mode equations 



3 d nin + 2) ^ t^/ \ 



/n(r,«:) = 



(5.12) 



and the following boundary conditions: 



fn (r, k) oc r 
fn{f, '^) oc exp(— /tr) / V Kr^ r 



oo 



(5.13) 



So /~ is regular at r = and is regular, i.e., bounded, as r ^ oo. 
For V{r) = these boundary conditions are those satisfied by In{r, n) and 
/C„(r, k), respectively. Furthermore, as the behaviour at r = is determined 
by the centrifugal barrier, and the behaviour for r — oo by the mass term, 
these boundary conditions are independent of the potential. If we write 



fn (r, k) = In{r, k) [1 + (r, k)] , 
f:{r,K) = IC^{r,^)[l + h+{r,^)], 



(5.14) 
(5.15) 
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then the functions h^{r, k) become constant as r — and as r — oo, and 
for finite r they interpolate smoothly between these asymptotic constants. If 
we impose, for r — >■ cxd the boundary conditions h^{r, k) —>■ the Wronskian 
of and /~ becomes identical to the one between Kninr) and /„(fi;r), i.e., 
equal to Applying the fluctuation operator to our ansatz, Eq. ()5.11|) . 

we then find 

1 1 °° 
+ ^2 ^ ^(^^ ^2) ^ _-5(r - r') ^(n + l)C^(cosx) 

n=0 

= \sir - r')-^6{x - x)^S{e - e')6{v - , (5.16) 
r-^ sm X sm U 

where we have used the addition theorem 

n I ^ 
E E ^n^-(^3)i;l^(f^3) = ^C^(COSX) (5.17) 

and the completeness relation for the 0(4) spherical harmonics. 
Numerically we proceed as follows: the functions satisfy 

which can be solved numerically. The second differential equation is solved 
starting at large r = f with /i^(f, k) = (f, k) = 0, and running backward. 
In principle we should take f = oo. However, iff is chosen far outside the 
range of the potential the fucntions h^{r, k) are already constant with high 
accuracy. In the numerical computation r = oo always means r = f with a 
suitable value of f. 

For the first differential equation we first obtain a solution /?.„(r, n) start- 
ing at r = 0, with hn{0, k) = h'^{0, k) = 0. This function does not satisfy the 
boundary condition required for the Green's function. The function h~{r, k) 
is obtained from hn{r,K) in the following way: from the definition of the 
functions h„ we have 



f-=Z,{l + h-). (5.20) 
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/~ and /~ are both solutions of the same hnear homogeneous differential 
equation and regular at r = 0, so they are proportional to each other, /~ = 
Cf~. The constant C follows from the boundary condition at r — ^ oo as 



C= . , (5.21) 

1 + K, (oo) 



and we obtain 



1 + hn{00, k) 

which obviously solves the differential equation with the appropriate bound- 
ary conditions. Of course in the numerical implementation r = cxd is taken 
as r = f (see above). 

Finally the Green's function is given by 



1 °° 
Ztt^ ^ — ' 

n=0 

(l + /i;(r<,z/^))(l + /i+(r>,i/2)). (5.23) 



In order to perform the subtractions needed in the process of renormal- 
ization, we not only need the functions which are exact to all orders in 
the potential V{r), but also the functions which are of first and second order, 
hn^"^ and hn^"^ , and the inclusive sums 

oo 

hF^ = E • (5.24) 

j=m 

This has been discussed in Refs. 1211 HI]- Obviously = hn^ as it 
includes all orders of V{r) except the zero order part. Writing the mode 
equations ()5.18p in the form 

Vh = V{l + h), (5.25) 

where T> denotes the differential operator on the left hand side, we have 

= V, (5.26) 
p/i(2) = , (5.27) 

Vh^ = Vh^ , (5.28) 
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and so on. Solving these and similar equations we find the parts of that 
correspond to a precise order in V{r). As V{r) is of a specific order in the 
couplings we thus may single out precisely specific perturbative parts. 

The rescaling needed in order to pass from h~ to h~ complicates affairs 
as it mixes orders. After some algebra one obtains 



and 



1 + hl^^' 



(5.29) 



(5.30) 



For the Green's function we likewise may define parts of a precise order in 
V{r); for x = x' we have 



g{x, x) = —^ ^(n + Iflnir, /€)/C„(r, «:) (5.31) 



X 



27T 

n=0 



1 + h^+'{r, k) + h^y{r, k) + h^^\r, K)h 



,(1) 



^(i)( 



1 



^(i)(x,a;) = V'(n + l)^X„(r, fi:)/C„(r, k 

n=0 



X 



,(1) 



^(2)(a;) = g(^)(a;,x) = —y2in + lYUr,^)ICn{r, 



n=0 



r, k) 



(5.32) 



(5.33) 



While Q^'^\x,x) is finite, Q^^\x,x) and Q^^\x,x) are divergent. Renormal- 
ization will be discussed in section |S] 



6 Computation of the Fluctuation Determi- 
nant 



The fluctuation determinant which appears in the rate formula 



V = det' 
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(6.1) 



can be written formally as an infinite product of eigenvalues of the fiuctuation 
operator. The prime denotes taking the absolute value and removing the 
translation mode. As in the previous section we introduce the generalization 

Note that we omit the prime, here. Using the decomposition of the Hilbert 
space into angular momentum subspaces we can write 

r ,.,2 _L ,,2 1 ^ r A/r 1 



n 




oo 

=n 






n=0 



with the radial fiuctuation operators 

(P 3d n(n + 2) 



r dr ^ 



det m'^\u^) 



+ + V{r) + 



(6.3) 



(6.4) 



dn denotes the degeneracy, in four dimensions we have dn = {n + 1)"^. 

According to a theorem on functional determinants of ordinary differ- 
ential operators [121 123 I2H] we can express the ratios of the partial wave 
determinants via 



det M„(z/2) 



^™ ,(0) 



detMi°V) ™^r(i^2^r) 
where ipni'^'^, r) and ilJn\i''^, r) are solutions to equations 

M„(z/2)V;.(z/^r) = , M.^^\u^)ij'^^\u\r) = 
with identical regular boundary conditions at r = 0. Of course 

V'f (z/^r)=X„(r,«:). 

Furthermore, we have 

V'„(z^^ r) = 1 + hn{r, k) J„(r, k) , 

where /i„(r, k) has been defined in the previous section. So we obtain 

det M^fz/^l 



det M^?V) 



1 + hJoo, k) 



(6.5) 



(6.6) 



(6.7) 



(6.8) 



(6.9) 
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and 



oo 



lnX>(z/^) = c?n In 1 + /i„(oo,k) 



(6.10) 



n=0 



The fluctuation determinant in the transition rate formula and in the 2PPI 
formahsm refers to the fluctuation operators at z/^ = 0, and so in the nu- 
merical computation we just need the functions /i„(0, oo), as for the Green's 
function. The only exception is the translation mode we will discuss in the 
next section. 

7 Unstable and translation modes 

In the one-loop formula for the transition rate the determinant of the fluc- 
tuation operator appears as det'(— (9^ + A^^), and the prime denotes two 
modifications with respect to the naive determinant: 

(i) the unstable mode has an imaginary frequency, corresponding to a 
negative eigenvalue ul = — z/^ of the fluctuation operator. It is to replaced 
by its absolute value. This mode appears in the s-wave n = and manifests 
itself by a negative value of 1 + /io(0, oo). So here we have to take the 
absolute value when computing the fluctuation determinant. There are no 
modifications of the fluctuation integral. 

(ii) the translation mode manifests itself, in the semiclassical approxima- 
tion, by the asymptotic limit 1 + /ii(z/^ = 0, oo) = 0. When backreaction is 
included, there is no exact zero mode, but a similar zero of 1 + hi{—uj^, oo) 
persists for a value of close to zero. We will identify it as the "would-be" 
translation mode, that would be again at cu^ = in the exact theory and, 
accordingly, we will treat it in analogy to the exact zero mode of the semi- 
classical approximation. Of course this is justified only as long as uf remains 
much smaller than the typical energy scales such as or u^. 

The fluctuation determinant ()6.2j) has a factor + u"^ which has to be 
removed, according to the definition of det'. Otherwise the logarithm of this 
expression, appearing in the functional determinant, does not exist. Further- 
more the Green's function is not defined either at u"^ = —w^. 

In the semiclassical approximation the translation mode is removed from 
the fluctuation determinant numerically in the following way JT]: we com- 
pute hi{oo, ±e^) for some sufficiently small e and replace in Eq. (j6.9|) 



1 + /ii(0, oo) 





2e2 



16 



i.e., we take the numerical derivative at uj^ = 0. 

In the backreaction computations we first determine the position of the 
eigenvalue by requiring l + hi{—ujf, oo) to vanish, and compute the numerical 
derivative not at i/^ = but at z/^ = —uj^, i.e., we remove a factor + . 

As the Green' s function is a functional derivative of the effective action 
we have to remove the zero mode from it as well. The Green's function in 
the n = \ channel has, at r = r', the form 

We can use the fact that the pole term is antisymmetric with respect to 
+ ijjI by computing the Green's function at = —oj^ ± and by taking 
the average of these two values. Then the pole term has disappeared and the 
averaged Green's function takes the form 

\ [^,(r, r, -uj + e^) + ^,(r, r, -u^j - e')] = RUt) , /^'^ ~ f , . 

n^O V In t ) 

(7.3) 

As long as uo"^ and are much smaller than the this is a good approxi- 
mation to the desired reduced Green's function 

[^i(-,-,0)U = 5^%^^. (7.4) 



In 



In the explicit numerical computation of the Green's function we use of 
course the expression 1)5.1111 . As evident from Eqs. ()5.22)) the pole arises 
from the rescaling of the mode function /ii(z/^,r), i.e., from dividing by 1 + 
/i_(i/^, oo). In averaging over the Green's functions at z/^ = —uj'j ± we add 
two very large terms which almost cancel. This can be done in a somewhat 
smoother way: if is sufficiently small we can assume that 1 + /ii(z/^,oo) 
passes through zero linearly and we may replace 

l + hi{-ujl ±e^,oo) ^ hi{-uj1 + e^.oo) -hi{-uol - e^,oo) . (7.5) 
The average over the Green's functions can then be cast into the form 

[yi\T-,^-,^)\red — r / 2,2 \ 1. I 2 2 N 

(7.6) 
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where /i(z/^,r) = Ji(/tr)[l + /ii(z/^,r)] is the mode function before the 
re-normahzation. 

As the translation mode is only approximate, the virial theorem men- 
tioned in section |21 no longer holds and we have to go back to the original 
expression for the prefactor, which is the normalization of the zero mode, so 
we compute the false vacuum decay rates via 

Of course the treatment of the approximate zero mode is an additional ap- 
proximation, beyond the one-loop or Hartree approximations. In the numer- 
ical computations the squared frequency of the zero mode is generally of 
the order 10"'^ in mass units. It becomes at most of order 10~^ near the 
critical points where our iteration ceases to converge, see section IHl So the 
approximation seems to be justified. 

8 Renormalization 

In the previous sections we have presented the basic formalism and its nu- 
merical implementation. There is still one point to be discussed: divergences 
and renormalization. As we are able to single out precise perturbative or- 
ders of the Green's functions and of the functional determinants, we can do 
the regularization of the leading orders analytically; we will use dimensional 
regularization. In particular we have 

and 

V = + ©(2) + p(3) ^ (8 2) 

where in both equations the first and second terms are divergent and the 
last term is convergent. The convergent parts are computed numerically 
using the methods described in the previous sections. The divergent terms 
are first computed and renormalized analytically. Then their finite parts 
are evaluated numerically, this is straightforward and involves some Fourier- 
Bessel transforms of the classsical profiles. 

Renormalization not only needs regularization but also renormalization 
conditions. In some previous publications jTUl dj we have used MS sub- 
traction which combines regularization and renormalization. Here we will 
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impose renormalization conditions which keep the effective potential close to 
the tree level one. For comparison we will also present results obtained in 
the MS scheme. 

The effective potential does not play any role in our computations, we 
always use the effective action. The effective potential is obtained if the 
classical field is homogeneous in space and time, or, as relevant here, in 
4-dimensional Euclidean space. The effective potential is often used when 
discussing quantum corrections, e.g., to Higgs potentials. It is relevant for 
the properties of the vacua, for double well structures it becomes complex 
around the maximum of the potential between the two wells. 

We now discuss the renormalization conditions. The false vacuum is 
relevant for the asymptotic region of the bounce as r — oo. It is convenient 
to require that in this region — > 0. This implies that the left minimum 
of the effective potential, for one-loop or Hartree back-reactions, remains at 
= 0, and that the mass remains the bare mass: 



It is reasonable, furthermore, to require that the vacuum expectation value 
in the true vacuum retains its tree level value, and that the energy difference 
e = — ?7(0tv) between the two vacua retains its tree level value, so that 
the semiclassical, one-loop and Hartree aproximations refer to essentially the 
same physical situation. Of course the exact shape of the effective potentials 
is different for these three cases. So we require 



The renormalization of the effective potential is discussed in Appendix 1X1 
For the one-loop back-reaction the counterterm potential is given by 



t/eff(0) 








(8.5) 



Using the definition 




7s ln47r - In — 

/i2 
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we find 



SP = + (8.7) 



2 



<5r/ = ISr/A^ + <5r^fi„ , (8.9) 
5A = -_L, + 5Afin. (8.10) 

Here the finite terms ^r/gn and 5Afin are the solutions of a linear system of 
equations given in Appendix ^ With these counterterms the equations of 
motion and the effective action of the bounce become finite. The relevant 
equations are given in Appendices IB. II and IB. 21 In the MS scheme the 
counter terms only consist of the parts proportional to and all the re- 
maining finite parts are set to zero. However, we have to deviate slightly 
from these conventions, 6p has to be chosen as in Eq. ()8.7p . otherwise the 
false vacuum is shifted away from = 0. The other finite parts have been 
set to zero. 

For the renormalization of the Hartree back-reaction it is essential that 
for a mass-independent regularization scheme all divergent parts are related 
[221 12n] . As a consequence these divergences can be conveniently removed by 
one counter term jSU] 

5C/di. = B{M^ - m') = -^(L, + 1){M^ - m^) , (8.11) 

see also the discussion in the Appendix of |31j. The part proportional to 
is an infinite renormalization of the vacuum energy. For the MS scheme this 
is all we have to do ^. If we want to impose the same boundary conditions 
as in the one-loop approximation with back-reaction, we have to introduce a 
set of finite renormalizations 

5C/fi„ = 5Afi, + 5pfi„0 + l6ml^4>^ - 5r/fi,03 + i^Afi^^^ . (8.12) 

/ o 



''in the strict sense B should be chosen equal to L^/16tt^ in the AIS scheme. This 
would lead to some tedious modifications of the back-reaction calculations. However, one 
can always change to if + 1 by modifying the renormalization scale fi. 
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Imposing again the conditions (j8.3|) and (|8.4p we find 



5Afin = 5pfin = W = . 



(8.13) 



Fixing the remaining counter terms becomes more involved; due to the non- 
hnearity of the gap equation we get a set of nonhnear equations. An iterative 
procedure is used to fix 5rj and 5\ numerically. This is discussed in Appendix 
IA.2I For the MS scheme the finite counter terms are set to zero, SUfiri = 0. 

With these counter terms the dynamical equations of the Hartree scheme 
become finite. The finite equations are given in detail in Appendix IB. 31 

9 Numerical results 
9.1 General remarks 

In discussing the numerical results it is convenient to use a parametrization 
which weights the relative importance of the classical and quantum parts of 
the action. One introduces the parameters 



potential, the mass m, we have set it equal to unity in our numerical com- 
putations. So all dimensionful parameters, like t] or e, and all dimensionful 
results are understood to be given in mass units. 

The parameter a parametrizes the shape of the potential, for small a the 
potential is strongly asymmetric, for a = 1 it becomes a symmetric double 
well potential. For values of a larger than 1 the role of true and false vacua 
is interchanged, so we can restrict ourselves to < a < 1. 

In the semiclassical approximation, i.e., one- loop without back- reaction 
one has 



dominated by the classical part and for small /3 by the quantum part. Of 
course this is not strictly true once one introduces back-reaction. Still we 
expect the effects of back-reaction to be small for large f3 and important for 



and the rescaling of the fields 



a 
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small (3. Furthermore, for large (3 the tunneling rate gets strongly suppressed. 
So our main interest will be in small values of j3. 

The self-consistent profiles 0(r) are obtained by iteration, each step of 
iteration consists in solving for a given J^{r) the equation for the profile (/>(r), 
and by computing then the new values of J^{x) for this profile. Of course 
in the first step J^{r) is not yet known; if the quantum corrections are large 
it is precarious to start the first iteration with J^{r) = 0. In order to avoid 
this problem we have computed, at fixed a, a series of solutions starting with 
a large value of (3, where the corrections are small. Then with descending 
values of [3 we have used in the first iteration step the self-consistent values 
of J^{r) of the previous value of (3. 

For the discussion of the data it is useful to introduce a shorthand notation 
for the four different cases we have investigated: 

- Case /: one-loop back-reaction, renormalization conditions Eqs. (j8.3|) 
and ()8.4|) . i.e. preserving the parameters of the true vacuum. 

- Case ir. one-loop back-reaction, MS renormalization. 

- Case nr. Hartree back-reaction, renormalization conditions Eqs. ()8.3|] 
and (Q 

- Case IV: Hartree back-reaction, MS renormalization. 

We display, in Fig. [T]the effective potential for a = 0.7 and (3 = 0.3 for 
the cases I to IV. The effective potential for cases / and /// is very close 
to the tree level potential over the whole range presented in the figure. The 
position of the true vacuum is seen to be shifted in different ways in the 
one-loop and Hartree cases. 

9.2 The bounce profiles 

In the semiclassical approximation the profiles 0(r) are, at fixed parameter 
a and arbitrary values of (3, determined by one universal profile 



which is independent of (3. Therefore, the change of the profile by the back- 
reaction can be displayed in a transparent way by plotting, at fixed a, the 
corresponding normalized profiles 



Mr) = P 




(9.4) 



0(r) = (3 




(9.5) 
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Figure 1: The effective potential for a = 0.7 and (3 = 0.3. solid line: tree 
level potential; dotted line: case /; dash-dotted line: case //; short-dashed 
line: case ///; long-dashed line: case IV] cases / and /// are hardly visible. 
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for different values of the parameter (3. For large /?, when the backreac- 
tion is weak, these profiles are expected to be independent of [3 and close 
to 0a(r). Indeed this is what we observe for [3 > 10. The /3-dependence 
observed for smaller values of 13 depends on the type of back-reaction and 
on the renormalization conditions. For the presentation in Fig. |21 we have 
chosen the Hartree approximation and the renormalization conditions which 
preserve the position of the true vacuum (case ///). With decreasing (3 the 
normalized profiles get lower and lower in the central region of the bounce. 
If (3 becomes smaller than 1, this decrease becomes substantial. For some 
lower limiting value of /3 ~ 0.2 the iterative procedure ceases to converge, 
during the iteration the profile collapses to 0(r) = 0. The qualitative be- 
havior is similar for all other cases and all parameters a we have considered 



In Fig. Olwe show the behavior of 0(0) / a//? near the critical values of 13 for 
the four different cases. This ratio becomes constant as /5 — > cxd and would 
be constant throughout in the semiclassical approximation. The behavior 
near the critical values is seen to be singular, as far as one can tell from a 
numerical computation. We therefore think that this phenomenon is genuine 
and not related to the way in which we iterate the equations. We have tried 
to modify the iteration by using the standard "overrelaxation" scheme 



where fi{r) is the i'th iteration of 0(r) and A^(r), and O is the operation 
that generates the new profiles in the normal iterative process (cr = 0). 
The runaway of the iteration below the critical value of (3 was persistent for 
various parameters a that we have tried out. The tendency towards = 
is apparent already for the parameter values of (3 where the iteration still 
converges, see Fig. |21 

For a < 0.6 and Hartree back- reaction the critical values of (3 are around 
(3 ~ 0.2 and then increase strongly. For the one-loop backreaction they 
increase roughly linearly from 0.2 to 0.6 between a = 0.3 and a = 0.7. At 
a = 0.8 the critical value is around /5 = 6 (!) for all cases. 

9.3 Effective actions and transition rates 

As mentioned above, the effective action in the semiclassical case is of the 
form 



(0.2, 0.3,... 0.8). 



/m(r) = ^/.(r) + (l-a)(0/,)(r) 



(9.6) 



'S'cfT — /5'S'cl[0] + 5'i_ioop , 



(9.7) 
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Figure 2: Behavior of 0(r) = 0(r)/^/^ for /3 = 20,1,0.5,0.25 and 0.2 at 
a — 0.6, for case III. In the central region, r < 4 these profiles are seen to de- 
crease with decreasing /? while they are independent of /3 in the semiclassical 
approximation. 
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Figure 3: Behavior of 0(0) / ^/J3 near the critical values which limit the region 
of convergence for a = 0.6. Short-dashed line: case J; full squares: case JJ; 
solid line: case /// and full circles: case IV. 
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where Sd and S'i_ioop are functions of a only. So we may display the effect 
of the back-reaction by plotting 



PSc\[(f)] + Si-loop 

where S^s is the effective action in the various approximations and renor- 
malization conventions. Si-ioop in the denominators depends on the renor- 
malization and is therefore different in the various cases. Clearly the ratio p 
should go to unity as /3 — >• cx), because there rj = 1/2a/5 and A = go to 
zero. Actually already for /? > 1 p is very close to one. Even near the point 
where the iteration ceases to converge, the deviation from unity is only a few 
percent. 

We display p in Figs. |3] and Figs. Both figures show the deviations 
for the four different cases. We see that the deviations for /3 > 1 are quite 
small; for (3 = 0.3 p is larger than 1 for all 4 cases, for (3 = 0.6 we have 
p < 1 for /5 > 0.5, for smaller values it increases strongly and displays 
some singlarity, a cusp or pole, near the critical value. The semiclassical 
action is dominated by the classical action, and so is the effective action 
including the various backreactions. We have, e.g., for a = 0.6 and case II, 
5'semi-ci = 395.08/3—29.522 where the first term is the classical action, and the 
second one the one-loop action. So the contribution of the quantum action 
is relatively small, and even substantial changes would not affect p. Indeed 
near the critical point it is the classical action which deviates strongly from 
its semiclassical value, and this is due to the strong changes in the profile 
0(r). We have to stress that here we are considering relative changes of the 
effective action. As the effective actions are of the order of a few hundred 
this implies absolute changes of several units, and this is what enters the 
transition rates. 

Another quantity of interest is, therefore, the ratio between the transition 
rates and their semiclassical values. In part this reflects the behavior of the 
effective actions, but also includes the ratio of the prefactors. For large f3 
these ratios will again go to unity, the logarithm will go to zero. We plot 
the logarithm of r/Fsomi-ci for small values of /3, for a = 0.3 and a = 0.6. 
We see that the ratio can amount to several orders in magnitude, for the 
one-loop backreaction. the rates are always suppressed with respect to the 
semiclassical rate. This is what we found as well in the 2-dimensional case JH] 
for the Hartree backreaction. Here we find, for the Hartree case, suppression 
for /5 < 0.5 and enhancement for larger values of a. 
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Figure 4: Ratio of Ses and 5*861111-01 as a function of P for a = 0.3; dotted line 
with full squares: case I; empty squares: case II; long-dashed line with full 
circles: case III; empty circles: case IV. 
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Figure 5: Ratio of Ses and 5*861111-01 as a function of (3 for a = 0.6; dotted line 
with full squares: case I; empty squares: case II; long-dashed line with full 
circles: case III; empty circles: case IV. 
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Figure 6: Logarithm of the ratio r/Fsemi-ci as a function of /5 for a = 0.3; full 
squares: case I; empty squares: case II; full circles: case III; empty circles: 
case IV. 
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Figure 7: Logarithm of the ratio r/Fsemi-ci as a function of /5 for a = 0.6; full 
squares: case I; empty squares: case II; full circles: case III; empty circles: 
case IV. 
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We have not been able to pinpoint precisely the reason for disappearance 
of the bounce solution at small j3. We can add the following comments: 

- the fluctuation integral shows no anomaly in the critical region. This 
means in particular that the instability of the classical solution is not 
caused by the appearance of a further unstable mode in the fluctuation 
operator. If the squared frequency of such a mode would cross zero, 
the fluctuation integral would receive a very large contribution propor- 
tional to the square of its wave function. We display the fluctuation 
integral a = 0.5 and for various values of (3 in Fig. |H1 In the semi- 
classical approximation it is independent of here we clearly see a 
/^-dependence which becomes stronger for small /3, but there is no sign 
of a singularity in this behaviour near the critical value (3 ~ 0.2. In the 
equation of motion for this integral is multiplied by A0 = so 
this contribution becomes more and more important for small (3. 

- the deviations of Ses from the semiclassical action are mostly due to the 
changes in the classical action, induced by the change of the profiles 
0(r). The one-loop and double-bubble contributions to the effective 
action do change substantially, but their influence on Ses remains nu- 
merically small. 

- the convergence of the partial wave summation both for the fluctua- 
tion determinant and for the fluctuation integral is very good; we have 
included terms up to n = 20 and included the asymptotic tail by an 
extrapolation (see Refs. PHICB])- The convergence remains excellent 
near the critical values of f3. 

- we have not spent too much effort in studying the region a > 0.7, where 
we approach the thin-wall regime. The classical actions become so large 
there that the semiclassical action is of the order of a few thousands 
(e.g. 1500/3 for a = 0.8) and correspondingly the transition rates are 
incredibly small. We note, nevertheless, that for a = 0.8 the critical 
values of (3 are around 6 in the various cases we have studied. 

We think that the instability of the bounce solution is mainly due to the 
term A0(r)jF(r) in the differential equation for the bounce. Its effects are 
difficult to analyze: we do not simply have a change of the effective potential, 
as displayed in Fig. [H but a modification that depends on r. 
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Figure 8: The fluctuation integral T{r) for the Hartree backreaction, case 
///, as a function of r for a = 0.5. The curves display a decrease with /3, 
which takes the values 10, 2.5, 0.9, 0.5, 0.3 and 0.2. 
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10 Summary and conclusions 



We have presented here two schemes for incorporating quantum backreac- 
tion into the computation of transition rates in quantum field theory, here 
apphed to false vacuum decay via a bounce: a one-loop backreaction, where 
the bounce profile is computed so as to minimize the effective action instead 
of the classical action; and the Hartree backreaction where the quantum back- 
reaction is included into the computation of the profile and of the quantum 
fluctuations. We have derived the general equations including their finite, 
renormalized form. We have used dimensional regularization and applied 
two different sets of renormalization conditions. 

We have presented numerical results for different parameter sets, using 
the a — f3 parametrization of Ref . [32] • The corrections to the semiclassical 
transition rate remain small as long as /3 = l/Arj'^ is greater than 0.6 and 
a = A/5 is smaller than 0.7. For lower values of (3 the deviations from the 
semiclassical transition rate become sizeable and can amount to several orders 
of magnitude. These deviations depend on the approximation, one-loop or 
Hartree, and on the renormalization conditions. The transition rates are 
reduced for the one-loop backreaction, for the Hartree backreaction they are 
reduced for a < 0.5 and enhanced for larger values of a. 

We find a critical value of (3, below which our iteration ceases to converge. 
The behavior of various quantities near this value suggests that this is a gen- 
uine phenomenon and not a technical deficiency of our iteration scheme. If 
this is so, it implies that bounce solutions do not exist below this critical value 
and that in this region the false vacuum decay proceeds via different config- 
urations in functional space. This is not unexpected for strong couplings: 
> 1 and A > a/0.3. 

Clearly, the computation of transition rates using the semiclassical ap- 
proximation and with quantum backreaction remain approximations which 
cannot be compared, at present, with experimental results or with lattice 
computations, such as those in Refs. [^3 EH IHS] • So it is not clear to what 
extent the inclusion of one-loop corrections, with or without backreaction, 
really represents an improvement. We may infer from our results that the 
semiclassical approximation is stable with respect to these higher corrections 
for a < 0.7 and (3 > 0.6, so we may believe it to be reliable. For lower values 
of f3 the deviations are strong, especially in the transition rates they amount 
to a few orders of magnitude. As the different approximations lead to quite 
different results, a comparison with lattice data would be of great interest. 
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The methods presented here can be apphed to various other computa- 
tions of quantum backreaction, e.g. to the selfconsistent computation of 
classical solutions which minimize the sum of classical and zero point ener- 
gies. While there one does not use the determinant theorem, one likewise 
employs techniques based on Green's functions in Euclidean space j^nHSHlEZl 
or techniques based on the analysis of phase shifts [SHlEnilini- Other possible 
applications of such selfconsistent computations include quantum corrections 
to instantons O E] , vortices jl2 113 El] and other classical solutions in 
quantum field theory. There are even cases without such a classical solution 
as, e.g., the selfconsistent pion cloud in the chiral quark model |lSllinillZlESl- 
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A Renormalization of the effective potential 
A.l The one-loop effective potential 

In the one-loop approximation the Green's function are computed with the 
effective mass 

= U"{4>) =m^- 6# + ^A02 . (A.l) 

Including the classical potential and the counterterms the 1-loop effective 
potential is given by 

f/cfT(0) = ^(0) + -lndet J' ^ +6Ui<P) 

= -^^ +8^ -64^(^-^^^ + 2j + 64^(^^ + 2 

2 8 

(A.2) 
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For the first and the second derivative of effective potential one finds 
U^g{(f)) = m2(/.-3#2 + ^(/)3 + 5p + (5m2(/)- 357702 (A.3) 

U'M) = -6#+-A(/)2 + 5m2- 65770 (A.4) 



(3A0-677)^ L,-ln-- -— — L,-ln-- + l 



We have already t/eff(0) — 0. The vanishing of the first derivative at = 
fixes 

The condition C/"(0) = ttt.^ leads to 

W = 5|,(iV + A™^)L, + ^. (A.6) 

The absolute minimum (true vacuum) of the classical potential occurs at 



If its position (the vacuum expectation value) is to be retained we have to 
require f/^jj(<;/)tv) = 0. If the energy difference e between the vacua is put 
equal to its tree level value we have to impose 

f/efr(0tv) = ^(0tv) . (A.8) 

The latter two conditions lead to 

= 97;A^ + 5r7fi,, (A.9) 
9A2 

5\ = ^L, + afin, (A. 10) 
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where the finite parts satisfy the hnear system of equations 
1 



1 



647r2 



3 Am 



2 i2 



(A.ll) 



(A.12) 



_ 6A0tv — 12?7 

with A12^ = 7W2(^^^)_ 



-•^tvan^-l)-m^ 



A. 2 The Hartree effective potential 

In the Hartree approximation the self-consistent effective potential as a func- 
tion of 4> is obtained from a variational potential 



u{(t)M^) 



t/(0) + -lndet ^ 



— A2 + 5;7div + 5[/fin(0) 



Afin + (^Pfin0 + 



-0 - (r^ + dr^fin)^ H 



In 



2 

3 



6477^ \ m 
+Adiv + SAt^ , 



m ( ^ 3 



8 

3A . 9 



(A.13) 



with 



by the condition 

The extremum is a maximum, so one has 



[/eff(0) =max U{(f>,M^) . 



(A. 14) 
(A.15) 

(A.16) 
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The condition ()A.15|) yields the gap equation 

3 

= m2 + 5m2^-6(r/ + 5r/fin)0 + -(A + 5Afin)0^ 

so that 



(A. 17) 



A^-^^L.-ln-^ + lJ+4i.A.^ (A.18) 

Nevertheless, for the variation of the potential L{{(f),A4'^) of Eq. ()A.13jl . A 
is defined as in Eq. flA.14|) . i.e., as a function of and Ai^. Cancellation of 
divergences in the gap equation requires B = L^/GAn"^ up to finite terms. It 
is convenient to choose 

L, + l 

Cancellation of divergences in the effective action then entails 



B = ^ . (A.19) 



Adiv = -Bm' = -^^m^ , (A.20) 

and we have 

A = f— In:^, (A.21) 

so that A(0,A^^(0)) = A(0,m^) = 0. If we want to impose, at </> = 0, the 
condition A^^(O) = we obtain 

5ml^ = . (A.22) 

The condition f/ofr(0) = fixes the finite cosmological constant to 

Afi„ = ^A2(0,m2) = 0. (A.23) 



We now consider the first derivative of the effective potential. It is given by 

with Ai^ taken as the solution of the gap equation. Thereby the second term 
vanishes and so 



UM) = <5pfin + m20-3(r/ + 5r/fi,)02 + ^±i^<^3 ^^.25) 
-I [Q{V + 5^fin) - 3(A + 5Afin)0] A(0, M^) . 
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The requirement C/gff(0) = 0, together with the condition A^^(O) = then 
leads to 

(^Pfin = 3(r/ + 5r/fi„) A(0, m') = . (A.26) 
With the present choice of the finite renormahzations we have 

(A.27) 

and 

Uls{(t>) = mV-3(r7+5r;fin)0'+^^|^0'-^ [6(77 + - 3(A + ^6^)0] A . 

(A.28) 

The condition f/cfr(0tv) = ^(0tv) yields 

(^%n0rv ^0tv = o n - - H ^ , (A.29) 

and the condition [/^^(^tv) = leads to 

?,5r^U<t>l + A) - ^0tv(0?v + 3A) = -i [677 - 3A0tv] A , (A.30) 

where in both equations and A are taken at = 0tv As ^T^fin and 
5Afin appear in the equations for M."^ and A we obtain a nonlinear system of 
equations for ^Tygn and (5Afin- We have solved this system numerically using 
an iterative procedure. 

In the MS scheme one generally lets — and docs not introduce any 
counterterms. Here we have to ensure that we have bare vacuum conditions 
at = with A1^(0) = w?. This is obtained by setting + 1 — > and 
omitting all further counter terms. Then of course, as in the one-loop case 
with MS scheme, the minimum at 0tv is shifted away from its bare value. 

B Renormalization of the effective action 

B.l Renormalization of the equation of motion 1-loop 
approximation 

The equation of motion for the classical profile is given by 

ct>" + -ct>' - U\4>) + 8U\ct>) + \u"\ct>)T = . (B.l) 
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The fluctuation integral is decomposed into the divergent leading order terms 
and a finite part as 

The computation of the finite part ^"^^^ has been described in section The 
leading order parts are given analytically as 

:F^'\x) = I T^T^ = (^e + 1) (B.3) 



and 



I i f d^kd^k' ^i(k-k')-{x-y) 
^1 (X) = - I d'yj (27r)8 (P + ^2)(^/2 + ^2)^(?/) 



(27r)4 y (27r)4(fc2 + m2)((A; + g)2 + m2) 
where we have defined the Fourier transformation 



The integration over k can be done and one obtains 



with 





2 + 4m2 









X 2 - In 



(B.4) 



V{q) = / d'ye-^'^yViy) . (B.5) 



•^^^n^) = -T^.L^Vix) + T^^ix) (B.6) 



a/I^P + 4m2 — |g| 

As the integrands, except for the exponentials, depend only on the absolute 
values of x and q the Fourier transforms reduce to Fourier-Bessel transform, 
we have 

47r2 1"^ 

V{q) V{\q\) = YT drr^Ji{\q\r)V{r) . (B.8) 
\^\ Jo 

and similarly for J-'^^\x) J-'^^\r). 
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In the equation of motion the fluctuation term and the counterterm po- 
tential are now given by 

SU'{(p) + l:U"'{(p)J^ = 5p + dm^(p - 35r](p'^ + ^^A^^ 



m 



;l, + i)-(-a02-6#) 



1 



167r2 



With the counterterms determined in Appendix A the divergent terms cancel 
and we get the flnite expression 



(B.9) 



+ l(3A0-6r;) 



B.2 Renormalization of the action in 1-loop approxi- 
mation 



The 1-loop part of the effective action is given by 

5i_/ = ^lnP. 



(B.IO) 



The logarithm of the fluctuation determinant can be expanded with respect 
to powers in the external potential V{x) as 



(-1) 



N 



(B.ll) 



-d' + U"{0) 

The flrst two terms in this expansion contain divergent parts. We write 

lnD = /l«--yl(2) + (lnD)(3y (B.12) 

and now consider the flrst two terms separatly. 

d^k 1 



Ad) 



(27r)4 F + m2 



d^xV{x) 



(B.13) 
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Using dimensional regularization we get 

= -^L, I dSV{x) + < (B.14) 



with 

For A^'^'^ we have 



2 



41 = / d'^xVix) . (B.15) 



d^k d^k' J d^xd^y e'^^-^'^-^'''y'^V{x)V{y) 
(27r)4 (27r)4 (P + m2)(fc'2 + m2) 

+ (B.16) 



with 



1 /""^ 



2 \/ (f ^ ^ a/q'^ + 4m2 + g 



9 \/ <f + — g 

"(B.17) 

With the counterterms determined in Appendix A the full one-loop action 
becomes 

•S'efr = Sc\ + Sii (B.18) 
= S^, + \\nV^+ j d'x (-5ryfi.03 + ^0^) - \Afl . 



The subtracted logarithm of the fluctuation determinant is evaluated accord- 
ing 
thei 
Eq. 



ing to sectional For y4g^\ as well as for J^^}, we have analytical expressions, 
their evaluation involves numerical Fourier-Bessel transforms like the one in 



B.3 Renormalization in the Hartree approximation 



Using the counterterms of Appendix IA.2I and the analysis of the fluctuation 
integral in Appendix IB. 11 the finite gap equation has the form : 



M\x) = m2-6(r/ + 5r/fi,)0(x) + -(A + 5Afin)0'(a;) 

+^A^fin(a;) (B.19) 
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with 

^fin(x) = + T'^'^Kx) + T^\x) . (B.20) 

The first term arises from JF^'') in Eq. ()B.3j) . the divergent part of JF^^) 
in Eq. ()B.6|1 and the counterterm ABM^ = M'^{L^ + l)/167r^ in the gap 
equation ()A.17|) . Once the profile and the fluctuation integral have been 
computed Eqns. ()B.19|) and ()B.20|) determine Ai^. The finite equation of 
motion becomes 

+ U'{<P) + 5UU(i>) - \{-Q{'n + %n + 3(A + 5Afi,0)^fi, = . (B.21) 
In this case the action is given by 

S'eff = Sel + 5Sfin + ^ In _ 1 ^(2) + ^ y d'^xV^x) 

- l>^fd'xJ^L (B.22) 



with 



SSfi^ = I d X 



-6v^^<P'{x) + ^-^<P\x) 



(B.23) 



In the MS scheme all finite renormalizations are omitted in the gap equa- 
tion, in the equation of motion and in the action, so that the latter reduces 
to 

Seff = Se, + i InV^ - + ^ I d^xV\x) 

- h [ d^xJ^l, . (B.24) 



Of course the fluctuations in InT'*^^^ and jFgn and the potential V{x) = 
Ai'^{x) — m?) are computed with a different M.^ and the profiles in A*^^) and 
in the classical action are computed using a different equation of motion. 
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